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Aims. The most important limitation for ground-based submillimetre (submm) astronomy is the broad-band absorption of the total 
water vapour in the atmosphere above an observation site, often expressed as the Precipitable Water Vapour (PWV). A long-term 
statistic on the PWV is thus mandatory to characterize the quality of an existing or potential site for observational submm- astronomy. 
In this study we present a three-year statistic (2008-2010) of the PWV for ground-based telescope sites all around the world and 
for stratospheric altitudes relevant for SOFIA (Stratospheric Observatory for Far-infrared astronomy). The submm- transmission is 
calculated for typical PWVs using an atmospheric model. 

Methods. We used data from lASI (Infrared Atmospheric Sounding Interferometer) on the Metop-A satellite to retrieve water 
vapour profiles for each site (11 in total, comprising Antarctica, Chile, Mauna Kea, Greenland, Tibet). The use of a single instru- 
ment to make the comparison provides unbiased data with a common calibration method. The profiles are integrated above the 
mountain/stratospheric altitude to get an estimation of the PWV. We then applied the atmospheric model MOLIERE (Microwave 
Observation and Line Estimation and REtrieval) to compute the corresponding atmospheric absorption for wavelengths between 150 
Jim and 3 mm. 

Results. We present the absolute PWV values for each site sorted by year and time percentage. The PWV corresponding to the first 
decile (10%) and the quartiles (25%, 50%, 75%) are calculated and transmission curves between 150 yt/m and 3 mm for these values 
are shown. The Antarctic and South- American sites present very good conditions for submillimetre astronomy. The 350 yt/m and 450 
Jim atmospheric windows are open all year long whereas the 200 //m atmospheric window opens reasonably for 25 % of the time 
in Antarctica and the extremely high-altitude sites in Chile. Potential interesting new facilities are Macon in Argentinia and Summit 
in Greenland that show similar conditions as for example Mauna Kea (Hawaii). For SOFIA, we present in more detail transmission 
curves for diff'erent altitudes (11 to 14 km), PWV values, and higher frequencies (up to 5 THz). Though the atmosphere at these alti- 
tude is generally very transparent, the absorption at very high frequencies becomes more important, partly caused by minor species. 
In conclusion, the method presented in this paper could identify sites on Earth with a great potential for submillimetre astronomy, and 
guide future site testing campaigns in situ. 
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1. Introduction 

Despite its large interest for astronomy, the submillimetre 
(submm) to far-infrared (FIR) wavelength range is mostly lim- 
ited by the pressure-broadened absorption of tropo spheric water 
vapou r present in the Earth's atmosphere (see, e.g. jMinier et alJ 
120101 for a summary of astronomy projects). The water vapour 
profile, and thus also the integrated column of water expressed 
in Precipitable Water Vapour (PWV), depends on the geographi- 
cal location. The PWV generally decreases with altitude, so that 
high-lying sites are best suited for submm- astronomy. Airborne 
astronomy, now possible on a regular base using SOFIA, is much 
less limited by water vapour and due to the lower atmospheric 
pressure, spectral broadening is less important. However, there is 
still absorption close to major and minor species (e.g., O3, N2O, 
CO) as well as collision induced non-resonant continuum ab- 
sorption. SOFIA can now access the very high frequency range 
(>2 THz), definitely closed for ground-based sites, but line in- 
tensities and abundances of the major and minor species need to 
be known precisely. First spectroscopic observations performed 



with SOFIA (see A&A special issue) showed that small-scale 
atmospheric variations in the tropospher e/stratosphere layer are 
very difl&cult to assess (iGuan et al ] l20T2h . 

For ground-bas e d sites, previous studies (e.g. 
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[ Schneider et al.l l2009l: fTremblin et al.l l201ll: [Matsushita etal] 
^9991: [Peterson et al.ll2003l) already showed that a few sites are 
well-suited for submm/mid-IR, and FIR-astronomy and their 
transmission properties are rather well determined (for example 
by Fourier Transform Spectrometer obervations in the 0.5-1.6 
THz range at Mauna Kea/Hawaii (iPardo et al.1 l200lh ). The 
high-altit ude (>5000 m) Chilean sites are known for d ry condi- 
tions (see [Matsushita et al.|[T999l:lPeterson et al.ll2003h . and site 
testi ng is now carried out at the driest place on Earth, Antarctica 
(see 'ChamberlmetaD 1 19971: I Yang etal.1 l201Ql: iTremblin et aP 
2011). Comparisons between Antarctic and Chilean sites are 
difl&cult and uncertain since they rely on ground-based instru- 
men ts that use differen t methods and calibration techniques 
(see lPeterson et al.ll2Q03l for example). The working conditions 
are also an important issue, a single instrument moved from 
Chile to Antarctica will have a different behavior in the harsh 
polar environment ( -70° C in winter at Dome C). A meaningful 
comparison is possible if several independent instruments are 
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used at each p l ace. A n example of such a study is the one of 
iTremblin etal](l201lh that obtained transmission data at Dome 
C thanks to radio-so undings and the radiometers HAMSTRAD 
(iRicaud et alJ l20ld) and SUMMIT08. However, it is rare to 
have many instruments at one site. The best solution is to use 
satellite data, which also enables to investigate any location 
on Earth. Thanks to the IAS I (Infrared Atmospheric Sounding 
Interferometer) on the Metop-A satellite, it is now possible 
to conduct such a comparison over several years with no 
instrument bias and with the same working conditions for the 
detectors. We present here a 3 -year study of the PWV of a 
selection of existing and upcoming submm- sites in Antarctica, 
Chile, Tibet, and Argentina, as well as for two SOFIA sta- 
tions, Palmdale/California and Christchurch/New Zealand. 
The individual and cumulated quartiles of PWV for each site 
allow a direct comparison. The transmission corresponding 
to these PWV values of the quartiles is then calcu lated using 
the atmospheric model MOLIERE-5 (lUrban et al.|[2b04) . The 
paper is organized as follows: After introducing all the sites 
that are considered in the present study (Sect. 2), we first show 
the PWV extraction from the satellite data and the comparison 
of the PWV statistics between all the sites (Sect. 3). Then we 
present the method to extract transmission curves with the 
MOLIERE atmospheric model and perform a site comparison 
for the 250 yum and 350 jim transmission (Sect. 4) The ratio of 
the monthly-averaged transmission to its fluctuations is then 
used to compare the stability of the transmission at 200 jim 
between all the sites at the end of Sect. 4. Section 5 deals with 
the results for SOFIA and Sect. 6 discusses and summarizes the 
results. 



2. Sites of interest 

We selected eleven representative sites around the world that 
can be grouped as follows: 

Antarctic sites 

Dome C: Concordia station with several telescopes, e.g. 
IRAIT (International Robotic Antarctic Infrared Telescope), 
'http://arena.oca.eu 

Dome A: new site in exploration, |http ://w w w. chinare . gov. cii 



South pole: SPT (South Pole Telescope)^ 

^http : //astro .uchicago . edu/research/south-pole-telescope| 

South- American sites 

Chajnantor Plateau in Chile: ALMA (Atacama Large Millimetre 

Array) site, |http://www.almaobservatory.org| 

Cerro Chajnantor in Chile: CCAT (Cornell Caltech Atacama 

Telescope), http : //www. ccato b servatory . org 

Cerro Macon in Argentinia 

Northern-hemisphere sites 

Mauna Kea in Hawaii: JCMT (James Clerk Maxwell Telescope), 

http : //w w w. j ach . hawaii . edu^ 

Summit in Greenland: http : //www.geosum mit.org| 

Yangbajing in Tibet: KOSMA (Cologne Observatory for 

submm-astronomy), h ttp : //www, astro .uni-koeln . de/ko sma 

"Stratospheric" sites 

SOFIA over Palmdale (US A) and Christchurch (New Zealand), 
|http : //www, sofia.usra.edu 



The location of the diff'erent sites are resumed in Fig.[T] Our 
selection indues known sites in Antarctica and North/South 
America plus northern-hemispheric sites that could be of 
interests (Summit in Greenland, Yangbajing in Tibet) to have a 
complete sky coverage. The geographical location, altitude, and 
related existing/planned telescopes are listed in Table [TJ 




South Pole Dome A 



Fig. 1. Location of the diff'erent sites of interest. 



Table 1. Location, altitude and telescopes for the diff'erent sites. 
Future telescopes are indicated in parenthesis. 



Site 


Latitude 


Longitude 


Alt. [m] 


Tel. 


Dome C 


75°06's 


123°23'e 


3233 


IRAIT 


Dome A 


80°22's 


77°2re 


4083 




South pole 


90° s 


0°e 


2800 


SPT 


Cerro Chajnantor 


22°59's 


67M5'w 


5612 


(CCAT) 


Cerro Macon 


24°3rs 


67°2rw 


5032 




Chaj. Plateau 


23°00's 


67°45'w 


5100 


ALMA 


Mauna Kea 


19°45'n 


155°27'w 


4207 


JCMT 


Summit (Greenland) 


72°35'n 


38°25'w 


3210 




Yangbajing (Tibet) 


30°05'n 


90°33'e 


4300 


KOSMA 


Palmdale 


34°37'n 


118°05'w 


12000 


SOFIA 


Christchurch 


43°31's 


172°38'e 


12000 


SOFIA 



3. PWV extraction from satellite data 

3.1. Method and results 

lASI (Infrared Atmospheric Sounding Interferometer) is an at- 
mospheric interferometer working in th e infrared, launched 



in 2006 on the METOP-A satellite dPhulpin et al.l 12007 



Pougatchev et al.l 120081: iHerbin et al.l 120091: IClerbaux et al.l 
2009|). The data are available at the website of the Centre for 



Atmospheric Chemistry Products and Service^ lASI was de- 
veloped by the French Space Agency CNESEl in collaboration 
with eumetsatEI. The satellite is on a polar orbit so that each 
point on Earth is seen at least once per day by the detectors. 
lASI is a Fourier Transform Spectrometer (FTS) working be- 
tween 3.7 and 15.5 jim. It is associated with an infrared im- 
ager, operating between 10.3 and 12.5 jim. Each pixel of the 
instrument corresponds to a spatial extent of 12 km at nadir, and 



^ http://ether.ipsl.jussieu.frj 
^ Centre National d' Etudes Spatiales 
^ European organization 
[http : //www ■ eumetsat . int i 



for meteorological satellites 



2 



Tremblin et aL: Worldwide site comparison 



vertical profiles of tropospheric humidity at ninety altitude lev- 
els (r esolution 1 km) are retrieved with a typically 10% accu- 
racy (iPougatchev et a l. 2008). The amount of precipitable wa- 
ter vapour is given by the integral of these vertical profiles. We 
averaged on a daily base all measurements whose central posi- 
tion falls in a zone of 1 10 km^ around each site to derive PWV 
statistics bet ween 2008 and 2010. This method was validated in 
lTremblin"era l. (201 1) in which we compared lASI satellite data 
with ground-based instruments (HAMSTRAD and SUMIT08) 
and radio- sounding data over the French-Italian base Concordia 
at Dome C in Antarctica. All instruments have a good correla- 
tion at low PWV which is precisely the range we are interested 
in for submillimetre astronomy. 



Dome C 




Time since 1st January [day] 




0.2 



^■ft.O 0.2 0.4 0.6 0.8 1.0 

Time percentage [%] 

Fig. 2. Top: PWV content measured at Dome C, Antarctica over 
the French Italian base Concordia between 2008 and 2010. The 
percentage indicates the fraction of days of the year on which 
we were able to extract the data from the satellite measurements. 
Bottom: function of repartition of the PWV for each year and for 
the whole period. 



Cerro Chajnantor 

6 1 . . . . 




Time since 1st January [day] 




ft.O 0.2 0.4 0.6 0.8 1.0 

Time percentage [%] 



Fig. 3. Same as Fig.[2]for Cerro Chajnantor in Chile. The satellite 
profiles were truncated at the altitude of the mountain summit. 
Note the diff'erent scale range of PWV compared to Fig.[2l 



The use of vertical satellite profiles is slightly trickier for a 
mountain site. Since we take all measurements in a zone of 1 10 
km^, we sometimes get profiles that do not contain mountain 
altitude but include lower ones. This would bias the retrieved 
PWV to high values. To overcome this difficulty, we generally 
truncated the profiles at the pressure le vel of the sit e of in ter- 
est. This method was already used by 'Ric aud et~aD (1201 Ol) to 
compare lASI measurements with the HAMSTRAD radiometer, 
over the Pyrenees mountains. They showed a very good corre- 
lation for the integrated PWV. We also determined in this way 
the PWV content at high altitudes (>11 km) over Palmdale, 
USA and Christchurch, New Zealand for the on-going and fu- 
ture flights of SOFIA. During a typical flight, SOFIA will range 
hundreds of kilometers from its base and the water vapour de- 
pends on local weather conditions at the position of the plane. 
Since the profiles are averaged on a 110-km^ area and the time 
statistics is done for three years, the results given here are repre- 
sentative conditions under which SOFIA is operating. Precise in- 
situ measurements of the P WV useable for ca libration are much 
more difficult to obtain (see lGuan et al 

The water vapour profiles are measured as a function of pres- 
sure by lASI. Therefore it is important to get the correct pres- 
sure level at the altitude of the sites. In- situ measurements were 
available for most of the sites to get the local averaged pressure 
level. Only for Yangbajing, Cerro Macon, and SOFIA, it was 
neccessary to use the pressure level of the US standard atmo- 
sphere (1976), computed at the altitude of the site. However, this 
method does not provide a precise pressure level just as a mea- 
sured one so that the retrieved PWVs are less reliable. By com- 
paring both methods for sites were measurements were avail- 
able, we estimate an error of around 25 % using the standard 
atmosphere. 

Figure [2] presents the PWV measurements for the Antarctic 
site of Dome C and Fig. [3] the measurements for the Chilean 
site of Cerro Chajnantor. The corresponding figures for all other 
sites are found in the Appendix and on a dedicated websit43- 
Since on Antarctic domes (expression for the mountains in 
Antarctica), the difference in height is small over 100 km, we 
do not need to truncate the profiles, whereas it has to be done 
for a mountain site like Cerro Chajnantor. Some ground mea- 
surements are available for different sites and can be compared 
to our results. Cerro Ch ajnantor and Chajnanto r Plateau were 
extensively studied (see iGiovanelli et al.l l200ll) and the PWV 
quartiles between 2006 and 2010 can be found in RadfordI 
(I2OIII) . The values deduced from opacity measurements at 225 
GHz are 0.61 mm (25 %), 1.08 mm (50 %), and 2.01 mm (75 
%) for the plateau; and 0.33 mm (25 %), 0.61 mm (50 %), and 
1.36 mm (75 %) for Cerro Chajnantor. The quartiles for Cerro 
Chajnantor are in a good agreement with our satellite measure- 
ments (between 10 and 20 % diff'erence). Note that the studied 
periods are not the same and can explain partly the difference. 
Furthermore, a scale height of the water vapour profiles on 
the Chajnantor area can be deduced from the quartiles of the 
plateau and the "Cerro". With the satellite measurements, we 
computed a 1/e scale height he of the order o f 1.1 km that is 
in a good agreement with measurements from IGiovanelli et aP 
(2001). The PWV quartiles at the South Pole were computed 
at saturati on from radio- soun ding balloons over a period of 40 
years (see IChamberlinll2002l) . These authors estimated that the 
PWV should be at 90 % saturated, thus the saturation level 
can be used as a proxy for the real PWV. They found averaged 
quartiles of 0.25 mm (25 %), 0.33 mm (50 %), and 0.44 mm 



ihttp://subirim.eu| 
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Table 2. First decile and quartiles of the PWV for all the studied sites. 



Time traction 2008-2010 


SOMA (Palm./Cnrist.) 


Dome A 


Dome C 


South Pole 


Cerro Chajnantor 


0.10 


f\ r\r\/^ ir\ r\r\ a 

0.006/0.004 


0.11 


0.17 


0.15 


0.27 


0.25 


0.006/0.005 


0.16 


0.22 


0.21 


0.37 


0.50 


0.007/0.006 


0.21 


0.28 


0.30 


0.61 


0.75 


0.009/0.007 


0.26 


0.39 


0.49 


1.11 


Time fraction 2008-2010 


Chajnantor Plateau 


Summit 


Cerro Macon 


Mauna Kea 


Yangbajing 


0.10 


0.39 


0.36 


0.47 


0.62 


1.21 


0.25 


0.53 


0.51 


0.66 


0.91 


2.47 


0.50 


0.86 


0.94 


1.02 


1.44 


inf 


0.75 


1.63 


1.96 


1.66 


2.57 


inf 



(75 %); the difference with our satellite measurements is around 
10 %. A radiometer at 225 GHz is operating at the CSCEl on 
Mauna Kea. PWV qu a rtiles b etween 1 9 97 an d 2001 can be 
found in lOtarola et al.l (l2010l) : iRadfordI (l201lh . 1.08 mm (25 
%), 1.82 mm (50 %), and 3.32 mm (75 %). The difference with 
the satellite measurements is around 15-20 % and can again 
be caused by variations between the period 1997-2001 and 
2008-2010. 

The effect of seasons are clearly visible in Figs. [2] and [3] 
Obviously, the PWV during winter is lower than in summer, 
i.e., below 0.3 mm at Dome C and below 0.6 mm at Cerro 
Chajnantor. The Chilean site presents a strong day-to-day varia- 
tion compared to the Antarctic site. This is caused by the daily 
variation of the PWV that is not present in Antarctica with the 
long polar night/winter. This may cause a slight bias in the 
non-polar site data because only one or two measurements per 
day are available, thus we do not sample very well these day- 
night variations. Nevertheless, we expect these statistics to rep- 
resent well the trends over years. The functions of repartition 
in Fig. [3] are very similar from 2008 and 2010, therefore the 
year-to-year statistics are not too much affected by the lack of 
the day-night variation sampling. In 2008 the PWV was slightly 
lower than in 2009-2010. This effect was seen in Antarctica 
(see Fig. [2]) also with the in-situ measurements (radio-soundings, 
HAMST RAD, SUMMIT08). P ossible explanations were dis- 
cussed in iTremblin etan(l201lh and should not be interpreted 
as a possible effect of the day-night variations. 

3.2. Site comparison 

Figure |4] shows the comparison between all sites as a function of 
repartition of the PWV over the period 2008-2010. Obviously, 
the air-borne SOFIA stratospheric instrument encounters the 
lowest PWV, lower than 0.01 mm most of the time and much 
lower than any ground-based sites. The Antarctic sites follow 
then as the driest places on Earth, with Dome A, and then Dome 
C and the South Pole. These sites have PWVs lower than 0.5 
mm most of the time (75 % of the time between 2008 and 2010). 
Chilean and Argentinian sites are next with a PWV lower than 
1.5 mm most of the time, followed by Hawaii and Summit in 
Greenland with a PWV lower than 2 mm, and finally Yangbajing 
in Tibet which presents a high level of water vapour. The satel- 
lite measurements saturate over the Tibet site with high values 
of the PWV (around 10 mm), this is clearly visible in the figure 
in the appendix. 

These results are presented quantitatively in Table[2lwith the first 
decile and the quartiles for each site. The difference caused by 
the truncation for mountain sites can be seen by comparing Cerro 

^ Caltech Submillimeter Observatory 



Cumulated : 2008-2010 




'•?).0 0.2 0.4 0.6 0.8 1.0 

Time percentage [%] 
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^•?).0 0.2 0.4 0.6 0.8 1.0 
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Fig. 4. Comparison of the functions of repartition of the PWV 
for all sites between 2008 and 2010. The percentage of days over 
this period on which we achieved the data extraction is indicated 
in parenthesis. The time percentage (x-axis) is made over the 
total number of days on which we achieved the data extraction. 
Note that in the upper panel, the curves for Palmdale (dark blue) 
and Christchurch (green) are blended and very close to the 0.0 
mm level. 



Chajnantor to Chajnantor Plateau. The mountain (Cerro) is close 
to the ALMA site, therefore the profiles that are extracted from 
the satellite are the same, the difference in the PWV retrieved 
is mostly caused by the altitude truncation of the profiles. The 
first decile and the quartiles again clearly show that Antarctic 
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sites are the driest sites followed by South- American sites and 
then northern-hemisphere sites. It is remarkable that our long- 
term satellite statistic of PWV shows that the site of Summit 
in Greenland offers comparable observing conditions (PWV and 
altitude) like the ones on Mauna Kea which opens a new per- 
spective for submm astronomy in the northern hemisphere. 

4. Atmospheric transmission from PWV and 
MOLIERE 

4.1. Method 




500 1000 1500 2000 

Frequency [GHz] 



3 mm 850 /j,m 450 f2.m 300 jj,in 200 150 jam 

Wavelength 

Fig. 5. Transmission curves between 150 fim and 3 mm corre- 
sponding to the first decile (grey) and the quartiles (red: 25 %, 
blue: 50 %, green: 75 %) of PWV between 2008 and 2010 at 
Dome C. 




500 1000 1500 2000 

Frequency [GHz] 



3 mm 850 /iin 450 /j,rn 300 /^m 200 iJ,m 150 jam 

Wavelength 

Fig. 6. Transmission curves between 150 jim and 3 mm cor- 
responding to the first decile (grey) and the quartiles (red: 25 
%, blue: 50 %) of PWV between 2008 and 2010 at Cerro 
Chajnantor. 



For the determination of the tropospheric transmission cor- 
responding to the PWVs of the various deciles and quartiles for 
each site, we use MOLIERE-5. 7 (Microwave Observation and 
Line Estimatio n and REtrieval), a forward and inversion atmo- 
spheric model (lUrban et al. 2004), developed for atmospheric 
science applications. It has previously been used to calculate 
the atmospheric transmission up to 2000 GHz (-150 yum) for 



a large number of astronomical sites. The results are published 
in [Schneider et al.1 (l2009h and on a websit^^. 

The code calculates the absorption of radiation in the mm- 
to far-IR wavelength range (equivalent from to 10 THz in 
frequency) considering the wet- and dry-air components of at- 
mospheric absorption. Spectroscopic lines such as those of at- 
mospheric water (H2O), oxygen (O2), and ozone (O3) absorb 
strongly at short wavelengths, while collisions of H2O with 
O2 and nitrogen (N2) result in continuous absorption across all 
wavelengths. For the line absorption, a radiative transfer model 
including refraction and absorption by major and minor atmo- 
spheric species is included (H2O, 02,63, NO2, HNO3, CO and 
other lines up to 10 THz). Spectroscopic parameters are taken 
from HITRAISQ and the JPL0 database. The wet and dry air con- 
tinuum absorption is calculated as described in [Schneider et al.l 
( 2009). Temperature and pressure profile s of the sites were deter- 
mined using the compilation described in lSchneider et al.l (l2009l) 
and given on the website. The colder Antarctic atmosphere pref- 
erentially populates the low-lying energy state of water, leading 
to a greater absorption for a given PWV than that of a non-polar 
site. This eff'ect is included in the model. 

We present here the transmission curves corresponding to the 
first decile and the quartiles of PWV for two sites of currently 
large interest. Dome C (Fig. [5]) and Cerro Chajnantor (Fig. O 
while transmission curves for the other sites are found in the ap- 
pendix. From Table [2] we extract that Dome C shows generally 
lower PWV values than Cerro Chajnator for longer time peri- 
ods. Comparing the transmission curves at the same PWV (for 
example 0.28/0.27 mm which corresponds to 50 % of the time 
at Dome C (blue curve) and 10 % of the time at Cerro Chajnator 
(grey line)) indicates, however, that the transmission is slighhtly 
higher at Cerro Chajnator. It corresponds to 20 (65) % at 200 jim 
(350 yum) for Cerro Chajnator with regard to -10 (60) % at 200 
yum (350 yum) for Dome C. The diff'erence becomes more im- 
portant for high frequencies (shorter wavelengths) because the 
absorption by water vapour an d the dry-air componen t are both 
a function of the altitude. See [Schneider et al.l (l2009l) for more 
details on the diff'erent absorption characteristics depending on 
altitude and wavelength. In any case, observations at 350 and 
450 yum are possible all year long for both observatories. 

4.2. Sites comparison 

Based upon the method described in iTremblin et al. we 
use MOLIERE to calculate the transmissions at 200 yt/m and 350 
Jim. We do not compute here the transmissions for SOFIA be- 
cause it is by far superior compared to all ground based sites. 
Its transmission curves will be discussed in a dedicated section 
(Sect. 5). 

Figure [7] shows the functions of repartition of the 350- and 
200- yum transmissions for all the other sites. In Table O we give 
the first decile and the quartiles of 200- and 350- yum trans- 
missions. 25 % of the time, a transmission better than 55-65 
% (10-20 %) at 350 yum (200 yum) can be found in Antarctica. 
For the first quartile, Cerro Chajnantor catches up the Antarctic 
sites thanks to its high altitude. However for the second and third 
quartiles Antarctic sites have a better transmissions based on the 
long term statistics, especially at 350 yum. For the three sites, 
Chajnantor Plateau, Cerro Macon and Mauna Kea, the transmis- 



^ http://subinin.eu 

^ High Resolution Transmission molecular absorption database 
^ Jet Propulsion Laboratory 
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Table 3. First decile and quartiles of the 350-yum (top) and 200-yum (bottom) transmissions for all the studied sites. 



Time fraction 


JJome c 


Dome A 


ooutn 


Cerro 


Chaj. 


Cerro 


Mauna 


Summit 


Yangbajing 


2008-2010 






Pole 


Chaj. 


Plat. 


Macon 


Kea 






0.10 


0.62 


0.72 


0.61 


0.65 


0.56 


0.49 


0.40 


0.42 


0.19 


0.25 


0.57 


0.67 


0.56 


0.58 


0.46 


0.41 


0.27 


0.31 


0.02 


0.50 


0.51 
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sion window at 350 jim opens significantly while it is only rarely 
possible to observe at 200 jim. 

The extreme stability of the Antarctic sites is clearly visible 
in Fig. [TJ There is a stable plateau for the transmission between 
the first and third quartiles that allows long time series and sur- 
veys to be conducted on these sites. This is the direct conse- 
quence of the stability of the PWV that was previously iden- 
tified at Dome C in Fig. [2j With its higher altitude, Dome A 
presents the best conditions on Earth for submillimetre astron- 
omy in terms of transmission and day-to-day stability. 

4.3. Transmission stability at 200 jim 

The stability of the transmission is an important parameter for 
large surveys and time-series studies. The PWV of Chilean Sites 
is highly variable from one day to another whereas Antarctic 
sites are much more stable (see Figs. [2l and [3]). The stability of 
the transmission leads to a high transmission even at long time 
fraction in the repartition functions of Fig. [71 In order to quan- 
tify and compare the stabili ty, we used the indicator introduced 
by iDe Gregori et al.l (l2012h . the site photometric quality ratio 
(SPQR). It consists in the ratio of the monthly averaged trans- 
mission to its monthly standard deviation, on a daily time-scale 



Cumulated : 2008-2010 



SPQR = {T)lcrT 



(1) 



The variations of this ratio for the transmission at 200 jim be- 
tween 2008 and 2010 are plotted in Fig. [8] and indicate that all 
temperate sites have a SPQR ratio lower than 1 while Antarctic 
sites have a ratio greater than 1 . Therefore, the fluctuations of the 
transmission are greater than the averaged transmission, on tem- 
perate sites. This quantifies the high variability of the transmis- 
sion at these sites. Note that the Arctic site on the Summit moun- 
tain is also highly variable, hence Antarctica is really unique 
even among polar environments. On the Antarctic plateau. Dome 
A has the best SPQR ratio with a monthly-averaged transmission 
that is typically 3-4 times higher than the fluctuations. Dome C 
achieves also very good conditions with a ratio of the order of 2- 
3 while the South Pole has a monthly-averaged transmission of 
the same order of the fluctuations comparable to the conditions 
reached at Cerro Chajnantor. 
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Fig. 7. Comparison of the functions of repartition of the 350- 
yum (top) and 200-yum (bottom) transmissions for all the sites be- 
tween 2008 and 2010. The x-axis time percentage is done over 
the period on which the data extraction for the sites was possible, 
the actual percentage of data we got is indicated for each site in 
parenthesis. 



5. Atmospheric transmission for SOFIA 

The first successfull science flights with SOFIA were carried out 
in 2011 and summarized in special issues of A&A (2012, Vol. 
542, spectroscopic observations with GREAT) and ApJ (2012, 
Vol. 749, continuum imaging). In particular the spectroscopic 
observations at high frequencies (e.g. the [CII] line at 1.9 THz) 



showed that although the transmission is very high at these fre- 
quencies (typical PWVs are in the range of a few [im up to 50 
yum for altitudes between 11 and 14 km), the calibration of the 
data using atmospheric models is not straightforward and re- 
quires a more detailed understanding of the atmospheric pro- 
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Fig. 8. SPQR ratio for the 200-yum transmission for all the 
sites between 2008 and 2010. The mean values between 2008 
and 2010 are Dome A 3.6, Dome C 2.7, South Pole 1.3, 
Cerro Chajnantor 1.1, Chajnantor Plateau 0.7, Cerro Macon 0.7, 
Mauna Kea 0.6, Summit (Greenland) 0.6, Yangbajing (Tibet) 
0.3. 



cesses at the tropo sphere/stratosphere border. As pointed out in 
iGuan et aP (l2012h , atmospheric models are particularly impor- 
tant because atmospheric absorption lines are more narrow at 
these altitudes and their s mall-s cale variation is very difficult 
to measure. In'Guan et al. [([2012h . three atmospheric models are 
compared: AM (iPaine.20 1 1|) . ATRA N (Lord 1992, priv. comm.) 
and MOLIERE (lUrban et aDl2004l) . Special emphasize is given 
separately to the wet- and dry-component. While all models 
agree well for the wet component, ATRAN does not include the 
collision-induced absorption by N2 and O2, and AM calculates 
a two times lower opacity than MOLIERE. Results of versatile 
models like MOLIERE depend much on choices made by the 
user. At this stage, it is not clear which models fits best for the 
data calibration and studies to compare the AM and MOLIERE 
codes are under way. 

In this paper, we present selected model results of 
MOLIERE calculated for 4 diff'erent altitudes (11 to 14 km in 
steps of 1 km), three PWVs (3, 10 and 20 yum), and a fre- 
quency range of 1800 to 2000 GHz. We chose this frequency 
range because it includes the astronomically very important 
atomic finestructure emission line of ionized carbon at 1900 
GHz. Figure [9] shows the transmission curves for the 4 altitudes 
at a PWV of 10 yum. Generally, the transmission is high (>80 
%) except of the deep and broad water absorption features and 
line absorption of minor species. Only in these ranges, the atmo- 
sphere is opaque even at those high altitudes. Figure [TOl shows a 
zoom into the frequency range around the [CII] 1.9 THz line for 
diff'erent PWVs. The transmission remains high (~75 %) even 
for a PWV of 20 yum. 




Fig. 9. Atmospheric transmission curve for 4 flight altitudes of 
SOFIA calculated with MOLIERE at a PWV of 10 yum. 



SOFIA 12km 




Fig. 10. Zoom into the transmission curve for the frequency 
range containing the [CII] 1.9 THz for 3 difl'erent PWV at an 
altitude of 12 km calculated with MOLIERE. The [CII] line is 
indicated by a dashed blue line. 



6. Summary and conclusions 

We presented PWV statistics and the corresponding transmis- 
sion curves for potential and existing submillimetre observa- 
tional sites all around the world. Thanks to satellite measure- 
ments (lASI) it is possible for the first time to conduct a compar- 
ison which is not biased by the diff'erent calibration techniques 
of ground-based instruments that have very diff'erent behaviors 
caused by their working conditions. This is particularly true for 
comparison between polar and non-polar sites because of the 
harsh environmental conditions at the poles. 

Comparisons were done for the PWV statistics and for the 
transmissions using the atmospheric model MOLIERE. For as- 
tronomical purposes, a comparison using the transmissions is 
mandatory since the altitude of the site has also a strong influ- 
ence on the resulting transmission. For a percentage of 25 % 
of the year, Cerro Chajnantor presents very good transmissions 
that are comparable to Antarctic sites. However on the long term, 
only Antarctic sites provide a stable transmission that will allow 
unique science cases to be studied there. Observations at 350 yum 
and 450 yum are possible all year long while the 200 yum window 
opens significantly for at least 25 % of the year. Furthermore, 
Antarctic sites are free of day-night variations thanks to the 
long polar night. All the other sites studied here have a highly- 
variable transmission. The typical monthly-averaged transmis- 
sion at 200 yum is lower than its fluctuations on all the sites ex- 
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cept for Dome A and Dome C. Therefore, these two sites will be 
unique for surveys and time-series studies in the submm range. 

The method used to compare the different sites is robust 
and based on only one instrument, lASI, and the atmospheric 
model MOLIERE. We derived in this paper the PWV and 
atmospheric transmission on well-known sites for submillimetre 
astronomy and showed that it is possible to retrieve statistics 
that are in a good agreement with in- situ measurements. A 
calculator to show these PWV statistics and to compute the 
corresponding transmission at any given wavelength is available 
to the communit}0 for all the sites presented here and for the 
three year 2008, 2009 and 2010. The year 201 1 will be added to 
the calculator. Other potential sites could be investigated upon 
reques{3- 
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In conclusion, this method could identify sites on Earth with 
a great potential for submillimetre astronomy, and guide future 
site testing campaigns in situ. 

Acknowledgements. We thank the anonymous referee for his/her comments that 
helped to improve this paper, in particular pointing out the importance of the 
pressure level to retrieve correct PWVs. 
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Appendix A: PWV statistics and transmission 
curves 




Frequency [GHz] 



Fig. A.l. PWV statistics (top) and transmission curves (bottom) 
for Dome A. The transmission curve for the first decile of PWV 
is in grey and the quartiles of PWV are given by red: 25 %, blue: 
50 %, green: 75 %. 
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Fig. A.2. PWV statistics (top) and transmission curves (bottom) 
for the South Pole. The transmission curve for the first decile of 
PWV is in grey and the quartiles of PWV are given by red: 25 
%, blue: 50 %, green: 75 %. 
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Fig. A.3. PWV statistics (top) and transmission curves (bottom) 
for Mauna Kea. The transmission curve for the first decile of 
PWV is in grey and the quartiles of PWV are given by red: 25 
%, blue: 50 %. 
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Fig.A.4. PWV statistics (top) and transmission curves (bot- 
tom) for Chajnantor Plateau. The transmission curve for the first 
decile of PWV is in grey and the quartiles of PWV are given by 
red: 25 %, blue: 50 %. 
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Fig. A.5. PWV statistics (top) and transmission curves (bottom) 
for Cerro Macon. The transmission curve for the first decile of 
PWV is in grey and the quartiles of PWV are given by red: 25 
%, blue: 50 %. 
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Fig. A.6. PWV statistics (top) and transmission curves (bottom) 
for Summit (Greenland). The transmission curve for the first 
decile of PWV is in grey and the quartiles of PWV are given 
by red: 25 %, blue: 50 %, green: 75 %. 
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Fig. A.9. PWV statistics for Christchurch New Zealand, at a 
stratospheric altitude of 12 km. 



Fig.A.7. PWV statistics (top) and transmission curves (bot- 
tom) for Yangbajing (Tibet). The transmission curve for the first 
decile of PWV is in grey and the quartiles of PWV are given by 
red: 25 %, blue: 50 %, green: 75 %. 




Fig.A.8. PWV statistics for Palmdale USA, at a stratospheric 
altitude of 12km. 
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